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Abstract. The potentiality of interpreting the D-meson decay-dynamics has revealed itself to be 
strongly dependent on our understanding of the light-meson sector The statistics collected by 
FOCUS is already at a level that manifests parametrization problems for scalar particles. In this 
paper the first application of the K-matrix approach in the charm sector is illustrated and preUminary 
results on the and Ds decays to three pions are shown. 



INTRODUCTION 

Charai-meson decay-dynamics has been studied extensively over the last decade. Dalitz 
plot analysis has indeed revealed itself a powerful tool for investigating the effects of 
resonant substructures, interference patterns and final-state interactions in the charm 
sector. The isobar formalism traditionally applied consists of a sum of relativistic Breit- 
Wigner resonances properly modulated by form factors (Blatt-Weisskopf barrier factors 
are generally used) and multiplied by an angular term assuring angular-momentum 
conservation. The potentiality to investigate charm dynamics has revealed itself to be 
strongly connected with a knowledge of the light-meson sector. In particular, the need to 
model the scalar particles populating our charm-meson Dalitz plots has led us to question 
the validity of the BW approximation for the description of a resonance. Resonances are 
associated with poles of the 5-matrix in the complex energy plane. It is the position of the 
pole that provides the fundamental, model-independent, process-independent resonance 
parameters. 

The fitting of data with simple Breit-Wigner formulae corresponds to the most el- 
ementary type of extrapolation from the physical region to an unphysical-sheet pole. 
In the case of a narrow, isolated resonance, there is a close connection between the 
position of the pole on the unphysical sheet and the peak we observe in experiments 
at real values of the energy. However, when a resonance is broad and overlaps with 
other resonances, this connection is lost. The BW parameters measured on the real axis 
(mass and width) can be connected to the pole-positions in the complex energy plane 
only through models of analytic continuation. A formalism for studying overlapping 
and many channel resonances was proposed long ago and is based on the ^-ma?nx|[lll2] 
parametrization. This formalism, originating in the context of two-body scattering, can 
be generalized to cover the case of resonance production in more complex reactions, 
with the assumptions that the two-body system in the final state is isolated and that the 



two particles do not simultaneously interact with the rest of the final state in the produc- 
tion process[3]. Its implementation allows us to include the positions of the poles in the 
complex plane directly in our analysis, embedding in our amplitudes the results from 
spectroscopy experiments flBl]. 



THE FIT FORMALISM 

The formalism traditionally applied to three-body charm decays relies on the so-called 
isobar model. A resonant amplitude for a quasi-two-body channel, of the type 

r + c 

^-^ a + b, 

is interpreted a la Feynman. For the decay D nnn of Fig. \^ a D ^ n current with 




FIGURE 1. The D+ KKK decay diagram. 

form factor Fo interacts with a di-pion current with form factor iv through an unstable 
propagator with an imaginary width contribution in the propagator mass. Each resonant 
decay function is thus, 

A=FDFr X \c\'\a\'Pj{co^@',,)^BW{mab) (1) 

i.e., the product of two vertex form factors (Blatt-Weisskopf momentum-dependent 
factors), a Legendre polynomial of order J representing the angular decay wave function, 
and a relativistic Breit-Wigner (BW). In this approach, already applied in the previous 
analyses of the same channels 10.13], the total amplitude (Eq.|2l) is assumed to consist 
of a constant term describing the direct non-resonant three-body decay and a sum of 
functions (Eq.[l]) representing intermediate two-body resonances. 

A(D) =floe'*+£«/^'^'A,-, (2) 

i 

Additional care has to be applied to describe the /o(980) resonance because of the 
opening of the KK channel near its pole mass. A Flatte coupled-channel parametrization 
Jalsh is generally adopted. 



The K-matrix formalism 



For a well-defined wave of specific isospin and spin IJ, characterized by narrow and 
isolated resonances the propagator is, as anticipated, of the simple BW form. In contrast, 
when the specific wave IJ is characterized by large and heavily overlapping resonances, 
just as the scalars, the propagation is no longer dominated by a single resonance, but is 
the result of complicated interplay among the various resonances. In this case, it can be 
demonstrated on very general grounds that the propagator may be written in the context 
of the K-matrix approach as 

{I-iK-p)-' (3) 

where K is the matrix for the scattering of particle a and b and p is the phase-space 
matrix. 

While the need for a K-matrix parametrization, or in general for a more accurate 
description than the isobar model, is questionable for the vector and tensor amplitudes, 
since the resonances are relatively narrow and well isolated, this parametrization is 
unavoidable for the correct treatment of the scalar amplitudes. Indeed the kk scalar 
resonances are large and overlap each other in such a way that it is impossible to single 
out the effect of any one of them on the real axis. 

In order to write down the propagator, we need the scattering matrix. To perform 
a meaningful fit to our three-pion data, we thus need a full description of the scalar 
resonances in the relevant energy range, updated to the most recent measurements in 
this sector. To our knowledge, the only self-consistent description of i'-wave isoscalar 
scattering is that given in the K-matrix representation by Anisovich and Sarantsev in 
llsj through a global fit of all the available scattering data from the nn threshold up to 
1900 MeV. 

To introduce the K-matrix formalism used in this analysis, we must recall some 
general, non-trivial, considerations. The production of an isoscalar, 5-wave state with 
an accompanying pion, involves, in the energy region relevant to this analysis, five 
channels, namely 1 = TtTt, 2 = KK, 3 = r]r] , A = r]r]' and 5 = multi-meson states (mainly 
four-pion states at ^/s < 1.6 GeV). The production amplitude in the particular channel 
(00++),-;r can be written as 

Fi = {I-iKp)r/Pj (4) 

where / is the identity matrix, K is the K-matrix describing the isoscalar s-wa\e scatter- 
ing process, p is the phase-space matrix for the five channels, and P is the 'initial' pro- 
duction vector into the five channels. In this picture, the production process is viewed as 
consisting of an initial preparation of several states, which then propagate via the term 
(/ — iKp)^^ into the final state. In particular, the three-pion final state can be fed by 
an initial formation of {KK)it, {KK)k, etc. The particular process we are dealing with, 
(;r+;r");r+, is then described by the amplitude F^. 
In this paper, the K-matrix used is that of ^ : 
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scatt 



and describe a smooth part of the K-matrix elements; the 



fee) 

The factor g\ is the coupling constant of the K-matrix ' pole a to meson channel 
the parameters ff^^ 

factor (;-y^Q)(i-^^Q) suppresses a false kinematical singularity in the physical region near 
the TZTZ threshold (Adler zero). The parameter values used in this analysis are listed in 
Tabled 

TABLE 1. K-matrix parameters. 
K-matrix 



mass poles 


Sn7i 


Skk 




8nn 


g-qr]' 


0.65100 


0.24844 


-0.52523 


0.00000 


-0.38878 


-0.36397 


1.20720 


0.91779 


0.55427 


0.00000 


0.38705 


0.29448 


1.56122 


0.37024 


0.23591 


0.62605 


0.18409 


0.18923 


1.21257 


0.34501 


0.39642 


0.97644 


0.19746 


0.00357 


1.81746 


0.15770 


-0.17915 


-0.90100 


-0.00931 


0.20689 


„scatt 

-3.30564 


/■scatt 
III 

0.26681 


/■scatt 

In 
0.16583 


/scatt 
/l3 

-0.19840 


/■scatt 
/l4 

0.32808 


/■scatt 
/l5 

0.31193 



1.0 



SAO 

-0.2 



These K-matrix values correspond to a T-matrix description of five poles whose 
masses and half-widths, in GeV, are (1.019, 0.038X (1.306,0.167), (1.470,0.960), 
(1.489,0.058) and (1.749,0.165); the analysis in B\ does not require the o. The 
decay amplitude for the D meson into three-pion final state, where TZ^TZ^ are in a 
(77^'' = 00^^)-wave, is thus written as 



-n 



X 



(a) 



'^m2 -m2 



prod lGeV2 



-F, = {I-iKp)-} 
-^T^ 1 ^ s-SA/2ml 

prod f {^s-Saq){1-Saq)' 



(6) 



where /3a is the coupling to the pole a in the 'initial' production process, /f™'^ and 
are the P -vector background parameters. In the end, the complete decay amplitude of 
the D meson into three-pion final state is: 



prod 



A{D) = aoe'* + £a,e'^'A/ + Fi 



(7) 



where the index i now runs only over the vector and tensor resonances, which can be 
safely treated as simple Breit-Wigner's (see Eq.[T]). In the fit to our data, the K-matrix 
parameters are fixed to the values of Table \\\ which consistently reproduce measured 
s-wave isoscalar scattering. The free parameters are those peculiar to the P-vector, i.e. 



Pa, ffj"'^ and ^o'"", and those in the remaining isobar part of the amplitude, at and di. 



prod 



Note that the K-matrix masses and widths do not need to be identical to those of the T -matrix poles in 
the complex energy plane. 
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FIGURE 2. D+ and Ds signals and sideband regions. 
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FIGURE 3. Three-pion D+ and D, Dalitz plots. 



THREE-PION PRELIMINARY RESULTS 

The three-pion selected samples (Fig.|2l) consist of 1527 ±51 and 1475 ± 50 events for 
the D+ and respectively. The Dalitz plot analyses are performed on yields within 
±2a of the fitted mass value. 



Ds — > 717171 



The Ds Dalitz plot (Fig.|3l) is clearly structured. Indeed, two /o(980) bands in the low 
and high projection at 1000 MeV, a band in the 1500 MeV high-projection region and an 




FIGURE 4. Ds Dalitz plot projections and isobar fit results 



accumulation of events in the comer of the folded plot are easily recognized. 

TABLE 2. Ds fit fractions and phases with 
the isobar model. 



resonance 


fit fraction(%) 


phase (deg) 


NR 


25.5 ± 4.6 


246.5 ± 4.7 


/2(1275) 


9.8 ± 1.3 


140.2 ± 9.2 


/o(980) 


94.4 ± 3.8 


O(fixed) 


5o(1475 


17.4 ± 3.1 


249.7 ± 6.4 


p''(1450) 


4.1 ± 1.0 


187.3 ± 15.3 



When the isobar model is applied, satisfactory fits can be obtained letting only the 
/o(980) and /o(1500) parameters float freely. We obtain mo = 975 ± lOMeV, T^k = 
90 ± 30MeV and Tkk = 35 ± 15 MeV for the /o(980) and mq = 1475 ± lOMeV and 
Tq = 112±24MeV for the /o(1500). The fit C.L., evaluated with a estimator over a 
Dalitz plot with bin size adaptively chosen according to the statistics is 11.5%. The fit 
results^ are reported in Table |2l 

All the established resonances decaying to k^k^ with a sizeable branching ratio are 
first considered in the fit. The set of contributions reported in Table prefers to fit coeffi- 
cients of more than 2a statistical significance. A fit with the /o(1500) parameters fixed 
to the PDG values has been also performed; the C.L. is 1% and the non-resonant fit frac- 
tion increases to about 40% (to be compared with the value in Table El), indicating how 
the flat non-resonant component, not having particular signatures, absorbs the residual 
parametrization problems. A large systematic effect should thus be associated with this 
channel decay fraction weakening the potentiality to establish a reliable estimate of the 
annihilation processes. The projection fit results are shown in Fig.|4l It is interesting to 



A fit fraction is conventionally defined as the ratio between the intensity for the single amplitude 
integrated over the Dalitz plot and that of the total amplitude with all the modes and interferences present. 



TABLE 3. Ds fit fractions and phases with 
the K-matrix model. 



resonance 


fit fraction(%) 


phase (deg) 


s — wave 


87.8 ± 1.7 


O(fixed) 


/2(1275) 


11.5 ± 1.2 


120.3 ± 5.4 


p"(1450) 


5.4 ± 1.1 


198.6 ±9.3 
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FIGURE 5. Ds Dalitz-plot projections with K-matrix fit results. 



note that the sum of the fit fractions is about 150%, pointing to a large, rather suspicious, 
interference effect among the contributions. 

One of the main advantages of the K-matrix approach is the straight-forward prescrip- 
tion for constraining the amplitudes to respect the two-body unitarity and the universality 
of the poles in the complex energy plane. The resonance parameters are fixed to the ref- 
erence paper |5i[]. The results are quoted in TableEI The contribution of the whole 5- wave 
in the Ds decay into three pions is computed in a single fit fraction. The fit C.L. is 3.2%. 
It is interesting to note here that the sum of the fit fractions is very close to 100%. The 
fractions of /2(1275) and p*^(1450) are quite consistent in the two approaches; the main 
difference concerns the scalar sector. A good fit to our Ds ^3tz signal can be obtained 
even without the non-resonant contribution. The K-matrix fit results are shown on the 
Dalitz projections in Fig.|5l 



The D+ — > TT^TT^n^ channel shows an excess of events at low tttt mass, (Fig. 
which cannot be explained with simple BW's corresponding to well-established reso- 
nances. A fit performed with a free parameter single-channel Breit-Wigner returns a 
mass m = 443 ± 27 and width F = 443 ± 80. The decay fraction and phases are reported 
in TablelH The /o(980) parameters are fixed to the values found in the Ds channel, where 




its signal is statistically more robust. The fit C.L. is about 10%; when the low-energy 
BW, o, is removed the C.L. drops to 10^^. 



TABLE 4. D+ fit fractions and phases with 
the isobar model. 



resonance 


fit fraction(%) 


phase (deg) 


NR 


9.8 ±4.3 


O(fixed) 


/2(1275) 


12.3 ±2.1 


-213.3 ± 17.7 


/o(980) 


6.7 ± 1.5 


-145.9 ± 17.7 


p(770) 


32.8 ±3.8 


62.9 ± 16.8 


5o(1475) 


1.8 ± 1.2 


242.3 ± 25.8 


/o(400) 


18.9 ±5.3 


- 96.9 ± 30.7 



When the K-matrix approach is applied to the D+, it is found that the P-vector degrees 
of freedom are able to reproduce the features of the low-energy nn distribution with a 
fit C.L. of about 3.5%. The K-matrix fit results on the Dalitz projections are shown 
in Fig. [71 the fit fractions and phases are reported in Table|5I 



TABLE 5. D+ fit fractions and phases with 
the K-matrix model. 



resonance 


fit fraction(%) 


phase (deg) 


s — wave 


66.5 ± 4.2 


101.8 ±22.5 


/2(1270) 


11.4 ± 1.4 


247.0 ±9.0 


p0(770) 


21.2 ±4.4 


O(fixed) 
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FIGURE 7. D+ Dalitz-plot projections with the K-matrix formalism. 

CONCLUSIONS 

The K-matrix approach has been applied for the first time to the three-pion Dalitz-plot 
analysis in FOCUS. The results are extremely encouraging since the same parametriza- 
tion of two-body resonances coming from light-quark experiments also works for charm 
decays; this result was not obvious beforehand. In the D+ izizn channel a good C.L. 
fit is obtained even without the o, suggesting that a lot of work has still to be done before 
solving the o puzzle. A global fit including the charm data and the analysis of the KKn 
channel could add important experimental information. The method explained here will 
find its full application to forthcoming excellent statistics of the charm experiments. 
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